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Abstract 

The proposed paradigm of plasmonic atoms and plasmonic molecules allows one to describe and 
predict the strongly localized plasmonic oscillations in the clusters of nanoparticles and some other 
nanostructures in uniform way. Strongly localized plasmonic molecules near the contacting surfaces 
might become the fundamental elements (by analogy with Lego bricks) for a construction of fully 
integrated opto-electronic nanodevices of any complexity and scale of integration. 
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I. INTRODUCTION 



A technological break-through in the fabrication of nano dimensional clusters and other 
metallic nanoparticles gave rise to a development of such nano-technological and nanoop- 
tical branch as nanoplasmonics, which is of great interest to physicists, chemists, material 
engineers, IT specialists, and biologists. The nanoplasmonics deals with conduction electron 
liquid oscillations in metallic nanostructures and nanoparticles, and an interaction of those 
oscillations with light (plasmon-polariton) [l, 2]. Such studies are also aimed at the oscil- 
lations of a crystal lattice (SiC, for example) in the nanoparticles, whose interaction with 
light has much in common with the plasmon oscillations (phonon-polariton) 3|]. 

The important feature of nanoplasmonic phenomena is the combination of a strong spatial 
localization and high-frequency (from ultraviolet to infrared) of electron oscillations. In its 
turn, strong localization leads to a giant enhancement of the local optical and electrical fields. 
These important features of plasmonic particles made it possible to discover quite a number 
of new effects. One of the most developed is the use of large local fields near plasmonic 
nanoparticles for enhancement of the Raman scattering cross-section. Recent experiments 
have shown that such an increase may achieve 10-14 orders of magnitude, which may help 
to resolve single molecules Qj-Q]. The local enhancement of the fields can also be used to 
increase the fluorescence intensity and to determine the structure of a single DNA strand 
without using the fluorescent labels Q-Q]. By using the nanoparticles of more complex shape 
one can provide enhancement of both the absorption and the emission of light by natural and 
artificial fluorophores lOj. On the other hand, the plasmon nanoparticles are proposed to 



be used in nanolasers 111 



and to stimu 



ate plasmonic oscillations in nanoparticles by means 
of the optical emission (SPASER) |l2|. Beside these new applications of the plasmonic 
nanoparticles, one can essentially increase the efficiency-cost ratio, for example, in solar 
batteries or light emitting diodes Q] by using the achievements in nanoplasmonics. 
And finally, it is awaited that the nanoplasmonics will make it possible to create a new 
element base (Genuine Integrated Optics and Ultracompact Optical Components) for the 
computers and data processing equipment by taking an advantage of small dimensions of 
metallic nanoparticles and fast speed of optical processes 15 1. 

An intricate spatial structure of the physical phenomena, which form the basis of 
nanoplasmonics, impedes a development of the latter. Very often, the numerical studies 
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do not allow one to explain the physics of the observed phenomena, while the analytical 
studies are mostly devoted to the case of spherical and spheroidal nanoparticles, which are 
very far from the synthesized nanostructures from the viewpoint of geometry and physics. 
Separate nanoparticles are often approximated by point dipoles. This permits one to ex- 
plain several phenomena (for example, a propagation of plasmons in a chain of nanoparticles 



16(|). But as a whole, the physics of plasmonic oscillations in complex nanostructures and 
nanoparticles remains insufficiently studied, and to make progress in this field one needs 
new approaches and new ideas. 

As one of such approaches we suggest that the plasmonic effects observed in clusters of 
nanoparticles should be described by the paradigm of plasmonic atoms and molecules, which 
has much in common with the paradigm of normal atoms and molecules. This would allow 
one to explain the complex phenomena of nanoplasmonics in terms of atomic and molecular 
physics, and to propose quite a number of new applications of the nanoplasmonics. 

Within the framework of our approach we propose to treat the plasmonic oscillations of 
the connected nanoparticles as the plasmonic atoms. The plasmonic atoms are analogous, 
in a certain way, to the electron oscillations at the nanocrystal quantum dots, which are also 



often called the artificial atoms 
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The bound states of plasmonic atoms in the clusters of nanoparticles are proposed to be 
called the plasmonic molecules, and finally, the plasmonic oscillations in the nanoparticle 
gratings are proposed to be called the plasmonic crystals. Such an approach seems to be 
efficient and useful as it follows from the strongly localized plas monic oscillations recently 



discovered in the clusters of spherical nanoparticles 18|, ll9|, |20| . These oscillations, being 
localized and of short-range, represent the plasmonic molecules, and on their basis it is pos- 
sible to create different kinds of plasmonic nanostructures and to develop new nanodevices 



(optical nanosensors of single molecules, nanoelectromechanical systems (NEMS 2JJ), and 
even plasmonic computers). 

II. PLASMONIC ATOMS 

Let us consider in more detail the analogy between the plasmonic and ordinary atoms. 
In the simplest case of a hydrogen atom, the energy levels are defined by the Bohr equation 
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E n = ~, n= 1,2,3,... (1) 

n l 

where _Ei=13.55 eV is the ionization energy of hydrogen atom. The wave function of an 
electron has the form, in this case: 



iP nlm (r,e,<p) = Z nl (r)Yr(6,<p), (2) 

where Z n i is the radial part of the wave function, and Y™, the spherical harmonic. For 
example, in the case of the first excited (2P) atomic state, the expression for the radial 
function takes the form: 



Z 21 (r) = N 21 I— J e , (3) 

where N 2 \ is the normalized constant, and r#, the radius of the first Bohr orbit. Function 
([3]) is represented in Fig. [TJ 

In the case of a spherically symmetric plasmonic atom, i.e. the plasmonic oscillations in 
a metallic spherical nanoparticle of the radius Ro, under the Drude dispersion law e (u) = 
l — (u)pi/uj) 2 (where u p i is bulk plasmon frequency), the spectrum of the plasmonic oscillations 
will take the form: 



/ n 

E n = huJ piy 2n + V n = 1 ' 2 ' 3 '--- ( 4 ) 
whereas the electric field potential, in respect to (jlj), will have the form corresponding to 
(T2]) with the radial wave function 

Z„ = N„ I " £ (5) 

\ (fl„/r)" +1 , r > Ro, 

where n — 1, 2, 3, ... and N n is some constant. Function ([5]) for the index n — 1 is illustrated 
in Fig. [T] by a dotted line. It is seen from this figure that wave functions of the ordinary and 
plasmonic atoms are similar. The spectra ([1]) and (jl]) also have much in common. 

The plasmonic oscillations have similar properties in the other spherically symmetrical 



layered systems (nanomatreshka) 22j and spheroidal nanoparticles (nanorice) [23]. 

The resonant plasmonic frequencies and potentials of plasmonic atoms in more complex 
and less symmetrical nanoparticles have a more complex form. For example, in the case of a 
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metallic nanoparticle having a form of a three-axial ellipsoid with the semi-axes a\ > CL2 > as, 
the plasmonic frequencies will be defined by the ratio 24] 



p K (01) g (01) 

n = 1, 2, 3, ... m = l, 2, ... 2n + 1, (6) 

where .E™ and F™ are the internal and external Lame functions j^], and a stroke denotes 
the derivative of function by its argument. The dependences of the resonant dielectric 
permeabilities for n = 3 are shown in Fig. [2] . The explicit expressions for eigen frequencies 
and eigen functions, e.g., for the case of n = 3 and m = 7, have the form (Guzatov, Klimov 
to be published) 



\ 



^37 — ^pl 

^37 = C 37 xyz, 
where C37 is a constant ; and I 



(aia 2 a 3 ) 



E 

v a=l 



! 123) 



(7) 
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(«+«i)(«+ a i)} 3A; 



The distribution of a surface charge for the plasmonic oscillations with n = 3 and m = 7 
derived from Eq. ([7]) has the form depicted in Fig. [3] (Guzatov, Klimov to be published). 



Even more complex properties characterize the cubic and polyhedra 



26(-[3l|. The spectrum of the cubic plasmonic atoms was studied in 



plasmonic atoms 



26(, and some of the 



first values take the form 



0.46225, 0.54473, 0.58722, 0.66372, 0.74953, 0.83918 



(8) 



Figure H] illustrates schematically the surface charge distribution for several states of a 
cubic plasmonic atom as found in [271 ] . The experiments studying the cubic plasmonic atoms 
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3l| proved to be in an agreement with the theoretical calculations. 



It is very important that both the plasmonic atom wavefunction (electric potential) and 
the field intensity are changing within the whole volume of a nanoparticle and the adjacent 
space, despite the fact that the structure of plasmonic atoms may be rather complex. By 
analogy with that, the surface charge in any plasmonic atom changes noticeably over the 
whole nanoparticle surface. Such a smeared spatial structure of plasmonic atoms is their 



characteristic feature, and can be used in quite a number of applications such as nanosensors 
of single molecules. However, as the distance between nanoparticles decreases they start to 
interact more intensively. And instead of the plasmonic atoms it is expedient to use the con- 
ception of plasmonic molecules to describe the plasmonic oscillations in the closely situated 
and strongly interacting nanoparticles. As in the case of ordinary atoms and molecules, the 
question about the existence of bound states of the plasmonic atoms emerges here first of 
all. 



III. PLASMONIC MOLECULES IN A TWO-SPHERE CLUSTER 



The evidence for the existence of bound states of the plasmonic atoms was first provided 



in 



18 



19| by the example of a cluster of two closely sitiated spherical nanoparticles (Fig. [5]). 



Generally speaking, the plasmonic oscillations in the cluster of two nanoparticles were 
studied earlier too {32]- {3?]]. However, it was assumed, both explicitly and implicitly, that 
the bound states are not formed in such a cluster, and the interaction assumes a scattering of 
plasmonic atoms on each other with a respective splitting of energy levels, and hybridization 
of wave functions (potentials) of single plasmonic atoms 0, That was the reason why 
the plasmonic molecules had not been predicted earlier. 

The discovered by us spectrum of the plasmonic oscillations in a cluster of two spherical 
nanoparticles is shown in Fig. [6] for the case of identical spheres (R\ = R2 = Ro). There are 
only symmetrical (T-modes) and antisymmetrical (L-modes) states of the scattering within 
the range of cu < uj p i/\/2 (e < —1). They exist at any distance between the particles, and 
continuously transit into the respective states of the weakly interacting plasmonic atoms 
with the characteristics (j4j) and ([5]), within large distances between the nanospheres. Just 
those scattering states had become the subject of investigations in [321]- [37] . 

In the region of u> p i > u> > tu p i/\/2 (0 > e > — 1) the plasmonic oscillations may exist 
only at small distances between the nanoparticles, and correspond to the bound states of 
the plasmonic atoms and molecules (M- modes). Note that the existence of the plasmonic 
molecules is by no means connected with a quasistatical approximation, and the calculations 
made within the framework of the total system of Maxwell's equations, with account for 
retardation effects, provide evidence for the existence of plasmonic molecules in the case of 
a small gap as compared with the radii of spheres (Guzatov, Klimov unpublished). 
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The analogous spectra appear for higher azimuthal numbers m, although there are some 
peculiarities in the spectra of the unbound plasmonic atoms. 

In the region of small distances between the nanospheres, the properties of the plasmonic 
molecules may be studied analytically within quasistatic approximation. For the Drude 
dispersion law, the spectrum of the plasmonic molecules takes the form: 



m u pi 

e£f = — (M + m + 8 m ) acosh (R 12 / {2R )) + ... 

M = 1, 2, 3, ... m = 0, 1, 2 ... (9) 

where R± 2 is the distance between the centers of spheres and 



5 = 1/2, S 1 = -0.08578, 5 2 = -0.2639, ... 5^ = -1/2, 

5 m = -1 - J- + J_ l — + (m > 0) (10) 

2 2m 8m 3 16m 5 128m 7 V ; V ; 

The potential of the axis-symmetric (m=0) state of the plasmonic molecule in a space 
between the spheres (— r] < rj < rjo) has the form: 



x/coshr/ - cos£ / 1 ^-^ cosh (On + 1/2) 77) 
1 " a (m g cosh ((»+ 1/2) , to ) P " (cOS ° 

cosh ((M+ 1/2)77) 



P M (cos , (11) 



cosh ((M+ 1/2) ?7o) 

where is the associated Legendre function and 77, £, the bispherical coordinates 



x 2 + y 2 + z 2 + a 2 ^ x 2 + y 2 + z 2 - a 2 

cothr/ = , cot 4 



2az ' 2a^x 2 + y 2 ' 

a = sjR 2 2 /A-Rl (12) 



and where coshr/ = R\ 2 / (2Po)- 

For the antisymmetric unbound plasmonic atom we have respectively 
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(L + m- 1/2)" 1 /acosh (R 12 / (2R )) + ... 



L 



1, 2, 3 



m = 0, 1, 2 ... 



(13) 



and 



"01/ 



X 



a/ cosh rj — cos £ 

sinh ((L + m — 1/2) rj) 

sinh ((L + m — 1/2) r/o) 



L+m-l 



tin 



(cos£) cos {rmp) . 



(14) 



In the case of large azimuthal numbers, m 1, one can also find simple asymptotic 
equations for the spectra of plasmonic oscillations in a cluster of two spheres. This case is 
of importance for calculation of van der Waals forces between nanoparticles. 

Consider in more detail a structure of the wave functions (electric potentials) of the 
plasmonic molecules. Figure [7] illustrates a wave function spatial distribution of a plasmonic 
molecule and unbound states of plasmonic atoms in the x-z plane. 

Spatial structure of antisymmetrical (L-modes) and symmetrical (T-modes) wave func- 
tions of the unbound plasmonic atoms, in the axis-symmetrical case (m = 0) corresponds 
to the structure of wave functions of the isolated atoms. Namely, a positive charge lies in 
one part of a sphere, while a negative charge, which is equal to the positive one due to 
the electrical neutrality of spheres, lies in the opposite part of the latter. An interaction 
between the plasmonic atoms is reduced, in this case, to a charge re-distribution on opposite 
semi-spheres. 

In the case of plasmonic molecules (symmetrical M-modes), the situation is the opposite, 
and both positive and negative charges are concentrated in a small region near a gap between 
the nanospheres. Far from the gap, the wave functions of plasmonic molecules practically 
vanish. 

In the case of m = 1, where the wave functions have an angular dependence cos (p or sin ip, 
the situation is analogous. The only difference is that the dipole moments of the unbound 
plasmonic atoms and molecules are directed along the x- or y-axis. 

As a whole, the wave function of the plasmonic atoms remains more or less distributed 
over the volume of both nanospheres, whereas the plasmonic molecules are strongly localized 
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near the gap between the nanospheres in a region with a characteristic dimension of the order 
of the gap size with the maxima at spherical surfaces. 

A strong localization of the plasmonic modes is also well seen in Fig. [HI which shows a 
distribution of the (surface) charge over a surface of spheres. To plot these distributions 
we have used a polar system of coordinates (p, 9) for each sphere in plane y = 0. The 
surface charge distribution was depicted with p = Rq (1 + a (9)) polar curve where a (9) is 
the dimensionless surface charge density as a function of polar angle of corresponding sphere. 

As the distance between the spheres increases, the localization of the plasmonic molecules 
decreases. At critical distance between the spheres the plasmonic molecules disappear, 
whereas the unbound atoms do not change. 

The localization of plasmonic molecules and unbound plasmonic atoms is different, which 
defines their different dependence on the exciting radiation. The polarizability of the un- 
bound plasmonic atoms is of the order of a nanosphere volume a ~ -Rq, and they effectively 
interact with the uniform external fields of the respective orientation and symmetry. And 
on the contrary, the plasmonic molecules have a comparatively small polarizability a ~ A 3 , 
where A is the value of the gap between the spheres. As a result, the plasmonic molecules 
are weakly excited by the uniform optical fields (as compared to the unbound plasmonic 
states). On the other hand, the plasmonic molecules interact effectively with strongly in- 
homogeneous fields that are localized near the gap between the spheres. Such fields are 
produced by the atomic and molecular radiation near the gap. This makes the plasmonic 
molecules extremely perspective from the viewpoint of a development of nanosensors and el- 
ements of the nanodevices, which are sensitive to single molecules. Quite probably, that the 
SERS from single molecules Q]- ?| is due to their interaction with the plasmonic molecules. 
It seems that hot points in fractal structures (Shalaev's hot spots are also due to the 
plasmonic molecules. 



IV. PLASMONIC MOLECULES IN MORE COMPLICATED SYSTEMS OF 
NANOPARTICLES 

Above we considered the plasmonic molecules and unbound states of the plasmonic atoms 
formed in a cluster of two identical nanospheres. However, the plasmonic molecules may 
be formed within a much wider set of clusters of the nanoparticles. The only necessary 
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condition for the formation of plasmonic molecules is smallness of the gap between the 
nanoparticles of the finite volume as compared to their characteristic dimensions. The 
plasmonic molecules may also be formed in connected particles with two or more nanobubbles 
or dielectric nanospheres inside them (Fig. [H]). 



A. Cluster made from different nanospheres 

First of all, the plasmonic molecules can exist in a cluster of two different nanospheres. 
Here the spectra and wave functions of unbound plasmonic atoms and plasmonic molecules 
are completely analogous to the case of identical nanospheres. If the distance between the 
spheres is small, then the plasmonic molecule spectrum may be written in the form: 



£i (u) - e 3 \ ( e 2 (w) -e 3 



El (U) + £3/ \B2 M + £3 

« exp [(2N + 2m - 1) fi] , 
N = 1, 2, 3, ... m > 1, (15) 

where the parameter Q can be found from the expression 

p2 _ p2 _ p2 

coshf]= 12 p- K \ (16) 

and where Ex, e 2 , £3 are the dielectric constants of the first and second spheres and the 
space between them, correspondingly. 

For the dispersion law of Drude for both spheres it can be found from f|T5|) 



2 (2 2 N 1 f I 2 2\2 

+^ ltl u 2 pl>2 exp [- (2N + 2m - 1) fi]} 1/2 , (17) 
where w p i t \ and tOpip are the bulk plasmon frequencies of spheres. 



B. Cluster made of two nonspherical particles 

Strong localization of plasmonic molecules allows one to suggest that their characteristics 
depend mainly on the radius of curvature of almost contacting surfaces and the distance 
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between the centers of the inscribed (or circumscribed) spheres. So, for any smooth non- 
spherical particle the properties of plasmonic molecules formed in the gap may be estimated 
by considering the plasmonic molecules in two spheres, which approximate the nonspherical 
particles in the contact point (see Fig. [TUT) . 

For the closely set semi-infinite bodies the situation turns to be more intricate, since in 
this case part of the charges may go arbitrary far from the contact area. Nevertheless, here 
exist a localized plasmonic oscillations, which are analogous to the antisymmetric unbound 
plasmonic atoms (L- modes) and plasmonic molecules (M-modes). 

In case of closely set semi-infinite hyperboloids (Fig. fTTa ) the variables in Laplace equa- 
tion can be separated within prolate spheroid coordinates [39|, and the localized plasmonic 
oscillation spectra have the form 



u m — Vpi {P'-i/2+iL m (^o) {P m \/2+iL ( _ ^o) ~ P-l/2+iL (Mo)) 

X (P-l/2+iL m (/^o) P-l/2+iL ( _ ^o) 

+P-l/2+iL m (~Mo) P-l/2+iL (Mo)) J j 
u m = u pl {P-l/2+iM m (Mo) (P-l/2+iM i~ + P-l/2+iM (Mo)) 

X (P-l/2+iM m (Mo) P-l/2+iM (~ A)) 

+ P-l/2+iM m ( — ^o) P-l/2+iM (Mo)) J , (18) 

where m — 0, 1, 2, ... and the indices L, M, which characterize the modes, continuously 
vary from to the infinity; /io = d/f (d is half of the smallest distance between the hyper- 
boloids and / is half of the distance between the hyperboloids' focuses). The corresponding 
wave functions (electric potentials) in the gap (— fi < \i < fi ) are 

# = P™ l/2+ iL (V) (P^L (/*) " P m i/2 +l L ("A*)) , 

iffi* 5 = ^-l/ 2+i M fo) (r i/2+!M M + P-l/2 + iM ("/*)) , (19) 

where the prolate spheroid coordinates 77, /1 are related with the Cartesian coordinates by 
the equations 
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(77-//)/ = + / 2 - 2 A, 

(77 + //)/ = y/R* + P + 2fz, 

Inside the hyperboloids the expressions for wave function of plasmonic atoms and plas- 
monic molecules are analogous to ( fl9l ). 

The axis-symmetrical (m = 0) spectra for several values of L, M =1, 3, 5, 10 are shown 
in Fig. [TTb as functions of distance between hyperbolic tips. The spatial distribution of 
corresponding wavefunctions (electric potentials) for L, M =1 are given in Fig. [121 

As seen from the figures, for both the symmetrical and antisymmetrical plasmonic os- 
cillations there is observed a strong localization of the modes, and this ensures an effective 
interaction with the ordinary atoms, molecules or quantum dots with plasmonic atoms and 
molecules near the gap between tips. One more interesting fact is that the spectra of these 
oscillations are quite close to the oscillations in a two-sphere cluster. This allows one to 
consider the localized plasmonic oscillations in the gap between semi-infinite bodies as the 
virtual plasmonic atoms and molecules. 



V. POSSIBLE APPLICATIONS OF PLASMONIC MOLECULES 



Obviously, strong localization of plasmonic molecules may be suitable for a number of 
applications. This is particularly important when one needs to ensure effective interaction 
between the nanolocalized light sources (molecules and nanocrystal quantum dots) and the 
nanoparticles and nanostructures, as well as the nano-electromechanical devices [2l|, in 
which the van der Waals forces play an important role 4(3, 41]. 



A. The effect of the plasmonic molecules on the van der Waals forces (V.V. 
Klimov, A. Lambrecht to be published) 

As is well known, the van der Waals forces are associated with spatial dependence of the 
vacuum fluctuation energy density. If the distance between the plasmonic nanoparticles is 
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small, then the main contribution to the van der Waals energy comes from zero (vacuum) 
oscillations of the plasmonic modes, i.e. the plasmonic atoms and molecules [42]. 




In case of two identical plasmonic nanospheres the contribution to the van der Waals 
energy comes from zero oscillations of the antisymmetric and symmetric unbound plasmonic 
atoms (L- and T-modes) and plasmonic molecules (M-modes) 



In spite of the fact that the expressions for the van der Waals energy of different modes 
are formally similar, the noted contributions result in totally different physical consequences. 
Figure [6] shows that the energy of antisymmetric unbound states of the plasmonic atoms 
grows with distance, and this results in the particle attraction. In a similar manner, the 
energy of the plasmonic molecules diminishes with distance, and thus results in the repulsion 
of the nanospheres. Symmetric states of unbound plasmonic atoms also lead to a very weak 
repulsion. 

The attraction or repulsion conditioned by different plasmonic oscillations is of a simple 
physical nature. First consider the axis-symmetrical oscillations (m = 0). Here the longitu- 
dinal antisymmetric modes are conditioned by the oscillations of the parallel dipoles which 
are attracting in accordance with the law of electrostatics (see Fig. [13] and Fig. [7j). In case of 
symmetrical oscillations the dipoles have opposite directions, and this results in the dipole 
repulsion (see Fig. [TBI and Fig. [7j). 

For m — 1, i.e. when the wave functions are described by the law of cos 93 or siny? (<£> 
is the angle describing the rotation around symmetry axis), the directions of the dipoles 
corresponding to different plasmonic states are changed. For antisymmetric modes the 
dipole moments in different spheres should be crosswise to the symmetry axis and should 
have opposite orientations (Fig. [13] and Fig. [7]). It is easily seen that this results again in 
the effect of attraction of the unbound plasmonic atoms. For the plasmonic molecules the 
dipole momenta should also be crosswise to the symmetry axis, but, due to the symmetry 




TjM 1 TjL 1 Tjl 

u vdW ' u vdW ' U v< 



vdW 
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with respect to z coordinate, should have the same direction (see Fig. [13]). This again leads 
to the repulsion of nanospheres with plasmonic molecules. 

Fig. O illustrates the dependence of the contributions from different plasmonic states 
into the van der Waals energy on the distance between the spheres. This was obtained 
by direct summation of all the plasmonic modes in (12T|) . As it was clear beforehand, the 
symmetric modes lead to the repulsion, while the longitudinal antisymmetric modes lead to 
the attraction of the nanoparticles. Unexpectedly, the repulsion contribution from plasmonic 
molecules is almost equal to the attraction contribution of the unbound plasmonic atoms. 
As a result, the total van der Waals energy grows with distance, i.e. it is the attraction 
energy, but the amount of this energy is by an order smaller than the energy to be obtained 
without taking into account the contribution from the plasmonic molecules. 

Direct measurement of the van der Waals force between the plasmonic nanoparticles may 
form the basis of the experimental proof of the plasmonic molecule existence. 

B. The effect of plasmonic molecules on the emission of ordinary molecules and 
quantum dots 

Since the plasmonic molecules particularly effectively interact with highly inhomogeneous 
electric fields, then one can develop different detectors of the ordinary single molecules or 
quantum dots. Fig. [15] illustrates the dependence of the radiation decay rate of spontaneous 
emission of different molecules, which fall within the gap between the two spheres, on the 
radiation wavelength. It is assumed that the nanospheres are made of Na with the plasmonic 
resonances in the visible range (Fig. fT5a ) or SiC with the phonon-polariton resonances in 
the IR range (Fig. fl5b) . The peaks in the right-hand part of the figures correspond to the 
interaction of an ordinary molecule and the unbound plasmonic atoms, while the peaks in 
the left-hand part correspond to the interaction with the plasmonic molecules. 

The analysis of the figures shows that, similar to the case of the van der Waals energy, the 
interaction with plasmonic molecules is of greater importance than that with the symmetric 
unbound plasmonic atoms. Of fundamental importance is the fact that in a homogeneous 
external field the plasmonic molecules are not effectively excited. This results in the follow- 
ing: the clusters of two and more nanoparticles (where the geometry allows one to excite 
the plasmonic molecules of a certain frequency), can be used in the development of effective 
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nanodetectors of single molecules of high signal/noise ratio, since the external fields at the 
plasmonic molecule frequency fail to excite the plasmonic molecule effectively. 

On the other hand, the effective interaction between the plasmonic and ordinary atoms 
and molecules can be applied both to the development of SPASER and single atom nanolaser, 
and the sets of such devices. 

C. Large cluster or infinite array of nanoparticles 

Probably, the most interesting consequence of the plasmonic molecule strong localization 
near the gap between the nanoparticles of arbitrary shapes lies in possible existence of plas- 
monic molecules in the clusters with more than two particles. Figure [16] illustrates a cluster 
of nine metallic nanoparticles, in which different configurations of plasmonic molecules may 
exist. In each of the twelve gaps a plasmonic molecule having the properties defined by the 
curvature radii of the closest surfaces and given by expressions (|9l)-(fTTT) may be excited. In 
this system one can also excite the pairs, triples, quadruples, etc., of plasmonic molecules. 
Under simultaneous excitation of several plasmonic molecules a weak interaction between 
the molecules may be easily taken into account using the perturbation theory. The account 
for such an interaction results in the hybridization of the wave functions of single plasmonic 
molecules and a corresponding splitting of the spectrum. Similar reasoning is also fully appli- 
cable to the large-size arrays of nanoparticles, and this allows one to speak about plasmonic 
molecular crystals. In other words, a strong spatial localization of plasmonic molecules and 
their weak interaction with each other allows an almost linear scaling of any simpler struc- 
ture of plasmonic molecules. To put it in simpler words, the plasmonic molecules allow one, 
similar to the Lego bricks, to produce the optoelectronic nanodevices of any complexity. 

In principle, the excitation of plasmonic molecules in the clusters and arrays may be per- 
formed with the help of homogeneous optical fields of a corresponding resonant frequency. 
However, such a method of excitation is not highly effective due to weak polarizability 
of plasmonic molecules. A positive feature of such inefficiency is a weak sensitivity of plas- 
monic molecules to the external noise fields. A more effective method to excite the plasmonic 
molecules lies in the excitation by means of spontaneous or stimulated radiation (SPASER 
[3])of usual molecules at a frequency of plasmonic molecule oscillation . In this case, de- 
pending on the orientation and location of the ordinary molecules ( or nanocrystal quantum 
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dots) with respect to the exciting field, the plasmonic molecules may be produced in the 
given part of a cluster or array, and, thus, to produce different initial conditions for the data 
processing and transfer. 

VI. CONCLUSION 

In this work the analysis has been made of the plasmonic oscillations in single nanoparti- 
cles and the clusters of nanoparticles. It has been shown that all the plasmonic oscillations 
may be classed by the degree of their localization. Such oscillations where the localization 
is mostly defined by the size of single nanoparticles we propose to name "the plasmonic 
atoms", and those strongly localized in the nanoparticle contact zones (without electrical 
contact) we propose to name "the plasmonic molecules", since they are the bound states 
of plasmonic atoms. It has been demonstrated that the plasmonic molecules almost do not 
interact between each other, and, so, can form the structures of arbitrary spatial complex- 
ity. The methods of plasmonic molecule excitation and its applications in the description of 
NEMS and the nanodetectors of single molecules are briefly discussed as well. 
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List of Figure Captions 



Figure [TJ Radial wave functions (a.u.) of usual and plasmonic atoms. The parameter 
R a = v2re or R a = R . The solid line corresponds to an ordinary atom and the dotted 
line, to the plasmonic atom. 

Figure [2]Plasmon frequencies for triaxial nanoellipsoid for n = 3 as a function of aspect 
ratio 03/01 for 02/01 = 0.6. Different colors denote different values m=l, 2, . . . 7. 

Figure [3] Surface charge distribution (a. u.) for a plasmonic atom in triaxial metallic 
nanoellipsoid for n = 3 and m = 7. The ratios of ellipsoid semi-axes 02/01 = 0.7 and 03/01 



Figure HI Normal modes (a.u.) corresponding to the six major absorption peaks of a 



Figure [5] Geometry of a two-sphere cluster. 

Figure [6] Spectrum of plasmon oscillations in a two-sphere cluster for m = 1 for different 
modes L, M, T = 1, 2, 3 as a function of distance between the spheres. 

Figure [7] Spatial distribution of electric potential of the modes with L, M, T = 1 in x-z 
plane, a, Case m — 0. b, Case m = 1. Radii of the spheres, Rq=50 nm; distance between 
their centers, _Ri2=105 nm. Surface of the nanospheres is denoted by a dotted line. 

Figure [8] Surface charge (a.u.) of different modes with L, M, T= 1, 2, 3 at m = 1. 
Radii of the spheres i?o=50 nm; distance between their centers, i?i2=105 nm. Surface of the 
nanospheres is denoted by the dotted line. 

Figure [9] Plasmonic molecules in the gap between nanobubbles in metal. 

Figure 1101 Approximation of nonspherical nanoparticles with equivalent spherical 
ones,which allows to estimate properties of plasmonic molecules. 

Figure 1111 Plasmonic molecules in the gap between two hyperboloids. a, Geometry of 
the problem, b, Spectra for m = and L, M — 1, 3, 5, 10. 

Figure 1121 Spatial distribution of wave functions for localized plasmonic oscillations in 
the gap between two hyperboloids for m = and L, M — 1. Distance between hyperboloids, 
2d=5 nm; the least curvature radius is 50 nm (cf. Fig. [7h). Solid line depicts surfaces of 
hyperboloids; dotted line corresponds to the cluster of two nanospheres with R = 50 nm 
and R12 = 105 nm. 



0.4. 
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Figure 1131 Illustration of attraction and repulsion of the nanospheres due to different 
plasmonic states. 

Figure 1141 Van der Waals energy of two-sphere cluster defined by different plasmonic 
states versus distance between the spheres. 

Figure 1151 Relative radiative decay rate of a molecule located between two nanospheres 
versus the radiation wavelength, a, Nanospheres made of sodium (Na). b, Nanospheres 
made of silica-carbone (SiC). In both cases, Rq = 50 nm and Ru = 105 nm. Direction of 
the molecular dipole moment is shown by arrow. The dielectric constant for Na is borrowed 

n n 

from [43j], for SiC, from [44 1. 

Figure 1161 Excitation of plasmonic molecules in clusters of nine nanospheres with usual 
molecules, a, Usual and plasmonic molecules are not excited, b, External field with hor- 
izontal polarization excite the usual molecules, c, Radiation of usual molecules excite the 
plasmonic molecules with horizontal polarization. 
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